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The effect of our Galaxy's motion through the Cosmic Microwave Background rest frame, which 
aberrates and Doppler shifts incoming photons measured by current CMB experiments, has been 
shown to produce mode-mixing in the multipole space temperature coefficients. However, multi- 
pole space determinations are subject to many difficulties, and a real-space analysis can provide a 
straightforward alternative. In this work we describe a numerical method for removing Lorentz- 
boost effects from real-space temperature maps. We show that to deboost a map so that one can 
accurately extract the temperature power spectrum requires calculating the boost kernel at a finer 
pixelization than one might naively expect. In idealized cases that allow for easy comparison to 
analytic results, we have confirmed that there is indeed mode mixing among the spherical harmonic 
coefficients of the temperature. We find that using a boost kernel calculated at Nside=8192 leads 
to a 1% bias in the binned boosted power spectrum at £ ~ 2000, while individual Ces exhibit ~5% 
fluctuations around the binned average. However, this bias is dominated by pixelization effects and 
not the aberration and Doppler shift of CMB photons that causes the fluctuations. Performing 
analysis on maps with galactic cuts does not induce any additional error in the boosted, binned 
power spectra over the full sky analysis. For multipoles that are free of resolution effects, there is no 
detectable deviation between the binned boosted and unboosted spectra. This result arises because 
the power spectrum is a slowly varying function of £ and does not show that, in general, Lorentz 
boosts can be neglected for other cosmological quantities such as polarization maps or higher-point 
functions. 



I. INTRODUCTION 

Forthcoming results from the Planck satellite will set 
new limits on the primordial CMB temperature field at 
small angular scales. It will also strengthen our knowl- 
edge about fluctuations at very large scales. In this 
regime, the temperature dipole is overwhelmingly the 
dominant contribution - a fact usually attributed to the 
peculiar motion of the solar system, the Galaxy and the 
Local Group through the universe. As pointed out in [I], 
the very high amplitude of this dipole (~ 100 times larger 
than the primordial fluctuations of all higher multipoles) 
induces a peculiar quadrupole at a level detected by the 
WMAP satellite and accessible to the Planck satellite. 
When extracted from its primordial counterpart, this in- 
formation can be used to further constrain the physics 
of our motion with respect to the CMB. We might also 
ask what are the contributions of this peculiar motion 
to the higher order multipole moments of the temper- 
ature power spectrum? This issue was first addressed 
in [2] , where it was shown that if the Cgs vary smoothly 
among neighboring £ modes in comparison to the aberra- 
tion kernel that boosts the CMB sky, then the main effect 
is a small and unobservable modulation of the temper- 
ature power spectrum. Furthermore, it was also noted 
in [2] that the multipolar convolution of the ai m s with 
the aberration kernel has an expansion which is in pow- 
ers of f3£, where j3 = v/c ~ 10~ 3 is the Lorentz boost 
parameter. From this it follows that the expansion in 
harmonic space is not guaranteed to converge rapidly for 
£ > 800, a region which most conspicuously concerns fu- 



ture measurements by Planck. Nevertheless, essentially 
all approaches devoted to deboost CMB maps were based 
on standard harmonic-space decompositions [5]-[5]. 

More recently, it has been proposed [Sj that one can 
deboost the CMB at the level of time ordered data. A 
follow-up paper showed that the power spectra derived 
from simulated aberrated maps affected cosmological pa- 
rameter estimation [7j. For their analysis they modified 
HEALPix code to produce a simulated, aberrated map 
from an input set of CgS. This analysis did not include 
a way to translate between aberrated maps and non- 
aberrated maps. Since CMB surveys directly measure 
time-ordered data which is converted to a temperature 
map, we will need a method for deboosting this temper- 
ature information for future analysis. In this work, we 
discuss our map deboosting method as well as resolution 
limitations and give an estimate for appropriate map res- 
olutions for high-^ CMB analysis. 

In this paper we assess the effect of a Lorentz boost on 
the temperature power spectrum through a full real-space 
analysis performed on temperature maps. There are two 
relativistic effects to consider: a Doppler-shift of CMB 
photons and an aberration. The Doppler shift of a map 
is determined the magnitude of the monopole and by the 
boost parameter. To lowest order in /3, the contribution 
is a dipolar function multiplying the whole CMB map. At 
small angular scales the corrections are dominated by the 
aberration effect, which translates, in real space, into a 
non-trivial distortion and overlap of CMB pixels. These 
effects in turn depend on both the resolution and the 
pixelization scheme used to cover the sphere. Given that 
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Lorentz boosts do not preserve pixel shape, pixels with 
equal areas but different shapes are distorted differently. 
Moreover, the overlap of boosted pixels with unboosted 
ones depends heavily on the tessellation of the sphere. 
Using the spherical projection pixelization scheme de- 
scribed in 8J, which maintains straight line boundaries 
between equatorial and polar cap pixels and is easily re- 
lated to the traditional HEALPix equal area pixelization 

we construct a real-space matrix that fully accounts 
for the pixel distortions and overlaps among pixels. Our 
approach not only circumvents complications associated 
with an expansion in f3 altogether, but also benefits from 
the already developed real-space tools used in standard 
CMB analysis. 

We organize the paper as follows: in Section [TT] we re- 
view the basics of boosting the CMB in harmonic space 
and stress its limitations. In Section III we present our 



general real-space approach and discuss its consequences 
when applied to the HEALPix pixelization scheme. We 
apply our method to simulated full sky maps and com- 



pare the resulting power spectrum in Section III we will 
describe the effect of map resolution on deboosting and 
present limitations in [Vj and we present concluding re- 
marks in section fVIIl 



II. 



LORENTZ BOOST EFFECTS ON CMB 
PHOTONS 



Relativistic effects have been outlined in [9], here we 
summarize relevant points and and describe our notation. 

If an observer in the rest frame of the CMB (denoted 
S) measures a photon of frequency v arriving along a line 
of sight n, then an observer in another frame, <S", that is 
moving with respect to the CMB at velocity vv will mea- 
sure the incoming photon to be arriving along a different 
line-of-sight, n', with a different frequency, v' . (Note 
that we will not concern ourselves here with any ambi- 
guities in determining S associated with the existence of 
inhomogeneities, in particular a cosmological dipole, on 
the assumption that that dipole is 0(1O -5 ) and smaller 
than the effects we will uncover.) The motion of the ob- 
server in S' thus induces two effects: a Doppler shift in 
the photon frequency and an aberration - a shift in the 
direction from which the photon arrives. These two ef- 
fects can be seen explicitly in the relation between n and 
h' 



cos 9 + P 
1 + P cos 6 



n — v cos a 
7 (1 + /3cos(9)' 



(1) 



where f3 = v/c and cos 9 = nv. The change in observed 
frequency in S" is given by a simple Lorentz transforma- 
tion 



v = 7^(1 + /3cos( 



(2) 



where 7 = (1 — /3 2 ) -1 / 2 is the standard Lorentz factor. 
This angle 9 is related to the angle 9' , measured in the 
frame 5", via 



cosfl' = 



cos 9 + p 
1 + P cos 9 



(3) 



For small /3 one can clearly expand cos 9' in a conver- 
gent series in /? around cos#. It would appear obvious 
that one can therefore do the same for the spherical har- 
monic function Yi m (9',(f>) oc P™(cos#'). However, this 
is not the case, the correct expansion parameter for Yi m 
is pi. This expansion will therefore break down (or at 
least converge poorly) when t > 1/P- (This becomes 
obvious when one recalls that Pi(cos9) has 0(€) zeros 
between —1 and +1, thus changing cos# by Oifi) moves 
one zero past the next for I > 2/p. ) This issue was orig- 
inally noticed in [2] where it was indeed recognized that 
this convergence breakdown leads to severe difficulties in 
evaluating the transformation of a spherical harmonic ex- 
pansion of the sky at I > 1/p. Nonetheless, they claimed 
that if the CVs are smooth functions of t compared to 
the transformation kernel, then this poor convergence is 
a red herring, and to second order in /3: 



C;«Q(l + 4/3 2 +G(/3 3 )). 



(4) 



We want to stress here that these assumptions are ide- 
alized, but far from ideal. In practice we do not measure 
the theoretical power spectrum, only its estimator (Ci). 
Real data is contaminated from many different sources 
which destroys the assumption of smoothness. This sug- 
gests a real-space approach is more appropriate since the 
aberration can be accounted for in terms of pixel distor- 
tions without the need to resort to a P£ expansion. We 
show in the next section how this can be implemented. 



III. THE REAL SPACE APPROACH TO BOOST 
CORRECTIONS 

We wish to transform the boosted, pixelized temper- 
ature map into an unboosted map that can be analyzed 
with techniques already available for CMB analysis. To 
do this we need to compute a boost matrix that charac- 
terizes how the unboosted pixels transform into boosted 
ones. 

The pixelized representation of the temperature fluc- 
tuation AT(n) can be defined by: 



AT(p) = ^- f AT(n)d 2 n, 



(5) 



where the integral is over the pixel p which has area J7 p . 
According to this definition, the value at pixel p' of the 
boosted fluctuation AT'(n') will be: 

ATV) = TT- I AT'(n')d 2 n' (6) 



http:/ /healpix.jpl. nasa.gov 



f D~ 1 (h)AT(n) d 2 h 
*'p' Jp' 



(7) 
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where we have made use of the following transforma- 
tions [9]: 

AT'(n') = D(n)AT(n) (8) 
d 2 n' = D- 2 {h)A 2 h (9) 
D(n) = 7(1 + /3v ■ n). (10) 

Note that when going from ffify to Q, the boundary of 
integration has not changed. In fact, the integral over p' 
should be replaced by a sum of integrals over all pixels 
overlapping with p' . Note also that the overlaps can (and 
do) occur among pixels that are not neighbors of the 
original pixel. We can therefore write 

AT'(p') = J2 A pp' D' 1 ( P )AT(p), (11) 
p 

where A pp i are matrix coefficients encoding the fraction 
overlapping of pixels. This transformation includes the 
pixelized versions of both Doppler, and aberra- 

tion, Appi, effects, and can be implemented entirely in 
pixel-space. The Doppler term is an overall dipole and 
significantly contributes only at very large scales. The 
boost matrix coefficients, Appi, on the other hand, are a 
function of pixels p and p' and strongly depend on the 
aberration effect. They contain the information about 
the fractional overlap of an unboosted pixel p with a 
boosted pixel p' . 

To compute the boost matrix, we make use of the 
spherical projection of the HEALPix sphere described in 
[5]. In this pixelization the pixel boundaries are straight 
lines in both the equatorial and polar-cap regions which 
makes the pixels easier to boost. To compute the ma- 
trix, we find which new pixel the original pixel center is 
boosted to and search its nearest neighbors and compute 
the overlap area of the boosted pixel using an approxi- 
mation of the exact area integral. These fractional area 
overlaps comprise the boost matrix. 

IV. COMPARISON WITH EXACT RESULTS 

To get a qualitative understanding of how a power 
spectrum is transformed under a Lorentz boost, we sim- 
ulated skies with power in only a single mode with 
i,m = 200, 0. We boosted this sky with /3 = 0.01, so 
that significant mode mixing would set in around I ~ 100, 
and extracted the resulting agj m i . The analytic calcula- 
tion of the integral for the boosted a' lm [TU] is tractable 
for this choice of multipole, and provides exact results for 
comparison to our numerical results. 

Figure [T] clearly shows that power is shifted from the 
mode which contained all of the power in the unboosted 
sky to nearby multipoles in the boosted sky. Addition- 
ally, we see that the analytic results are in good agree- 
ment with the numerical values for an Nside=512 map. 
The discrepancies between the numerical and analytic 
calculations can be attributed to errors induced by pix- 
elization of the temperature map. We will discuss this 
effect at length in the following section. 




FIG. 1. Comparison of boosted (blue dashed line) and un- 
boosted (black solid line) temperature power spectra from a 
simulated sky at Nside=512 with all power in the I = 200, 
m = mode to the exact solution. Differences between the 
numerically boosted result and the exact solution (magenta 
dots) are predominantly due to map resolution effects. A /? 
value of .01 was used to obtain these results. 



V. RESOLUTION EFFECTS 

We have found that the resolution of the input map 
used for deboosting has a direct effect on the resultant 
power spectrum. This is due to pixelization effects - 
the inaccuracies inherent in replacing the true mapping 
between points on the sphere that the boost represents, 
with mappings from pixels of the unboosted sky, to pix- 
els of the boosted sky. The fractional overlaps contribut- 
ing to a pixel average inevitably includes area that the 
original, boosted map should not include. It is therefore 
important to check that deboosting is done at a suitable 
resolution for whichever range of t is to be included in 
the subsequent analysis. Higher multipoles are most sus- 
ceptible to pixelization-induced effects. If the map reso- 
lution is sufficient, then after a map has been deboosted, 
an equal boost in the opposite direction should repro- 
duce the original power spectrum over the full range of 
I to the desired/required accuracy. This can be used to 
check for residual resolution effects in the map at any 
preferred scale. 

Figure [2] shows how resolution effects change the power 
spectrum for maps at resolutions of Nsidc=2048, 4096, 
and 8192, which are higher than those traditionally used 
for CMB analysis. We show that in order to circum- 
vent the problem one must overpixelate the boost kernel 
for accurate deboosting. The Nside=8192 map, which 
has been deboosted with a boost kernel calculated at 
Nside=8192, has a 1% difference between the unboosted 
and boosted power spectrum for I ~ 2000. We confirm 
that this is a resolution based effect, since an Nside=8192 
map which has been boosted then deboosted has a 2% 
residual difference in the power spectrum for I ~ 2000 
as seen in Figure [3j in agreement with the conclusion 
that the 1% inaccuracy in the power spectrum in Fig. [2] 
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FIG. 3. Comparison of power spectra produced from a map 
that was not boosted versus a map that had been boosted 
and deboosted. The difference should be zero, as deboosting 
should return the map pixels to their original values. The 
difference shown is therefore due to pixelization effects. The 
blue line shows the difference for individual Ces and the black 
line shows the difference in binned Ces. Values are quoted in 
percent. 





FIG. 4. Percent difference between a power spectrum gen- 
erated by CAMB [TT] with n s - la and a standard LCDM 
power spectrum. Values are quoted in percent. 



FIG. 2. Comparison of power spectra produced from 
boosted versus unboosted full-sky maps at Nside=2048(top) 
4096(middle) and 8192(bottom). The blue line shows the dif- 
ference for individual Ces and the black line shows the differ- 
ence in binned Ces. Values are quoted in percent. 



is predominantly due to resolution effects. Interestingly, 
the fluctuations in the individual C#s are small for the 
deboosted map, suggesting that the ~ 5% fluctuations of 
boosted CgS around the binned bias in Fig. [2] are not due 
to pixelization effects. 

These resolution effects can clearly lead to inaccuracies 
in the extracted cosmological parameters if the difference 
in produced power spectra are assumed to be primordial 
in origin rather than due to the data processing. 

We show in Figure[4]the percent difference between two 



power spectra that have different values of n s (one was 
produced with the WMAP 7- year best-fit value, the other 
power spectrum was produced with the best-fit value plus 
one sigma using CAMB^dU ), showing that a difference 
in the spectrum at large scales gives a similar difference 
as the resolution effects at Nside=4096. This suggests 
that resolution effects could be confused with changes in 
cosmological parameters unless one is careful to deboost 
with the appropriate boost kernel. As shown in Fig. [2] 
an Nside=8192 boost kernel would be sufficient to see 
changes in parameters that result from deviations larger 
than 1% in the binned power spectrum. 



http: / / camb.info 
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FIG. 5. Comparison of Nside=2048(top) and 

Nside=4096(bottom). These have been boosted and 
downgraded to Nside=512. Values are quoted in percent. 

A recent paper by Catena and Notari [7J claimed that 
n s extracted from simulated maps at a resolution of 
Nside=2048 with intrinsic aberration and WMAP-like 
noise differed by approximately 1 sigma from value ex- 
tracted from an unaberrated map. 

We also note that boosting at a high resolution and 
downgrading the map before extracting the power spec- 
trum gives better results than simply boosting at the 
desired resolution from the onset as seen in Figures [5] 
& [U Therefore one should calculate their boost kernel 
at a pixelization well above the map pixelization before 
applying it to a map. 



VI. ANALYSIS USING MASKED MAPS 

Because masking a map is known to correlate modes 
at all scales [T2], we wanted to investigate the effect of 
both boosting and masking a temperature map on the ex- 
tracted power spectrum. To that end, we employed two 
methods. We used boosted maps at Nside of 2048, 4096, 
and 8192, rotated them to the direction of the dipole, 
downgraded them to both Nside=512 and 1024 (resolu- 
tions traditionally used for CMB analysis), and masked 
then using the conservative galactic mask released by 
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FIG. 6. Comparison of Nside=2048(top) and 

Nside=4096(bottom). These have been boosted and 
downgraded to Nside=1024. Values are quoted in percent. 

WMAP from the LAMBDA website El The same was 
done for unboosted maps for a comparison between the 
resultant power spectra shown in Figure [7] for Nside=512 
and Figure [8] for Nside=1024. We confirmed before do- 
ing full mask analysis that rotating each map did not 
contribute to our results. 

We also directly masked Nside 2048 and 4096 maps 
with a 20 degree wide isolatitude galactic cut. This was 
done because the released galactic masks are not pro- 
duced at resolutions higher than 1024. The results are 
shown in Figure [9j 

For the high resolution maps with isolatitude masks, 
we have found that masking a map produces a negligi- 
ble effect on the power spectrum extracted from a boost 
map. The most noticeable difference is that individual 
CiS oscillate by a wider range than the unmasked power 
spectra, but the binned differences between masked and 
unmasked are nearly identical. 

For the maps that have been masked with a galactic 
cut, we see that boosting at a much higher Nside, down- 
grading the map and then masking at either Nside=512 



s http://lambda.gsfc.nasa.gov 
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FIG. 7. Comparison of nside=2048(top), 4096(middle), and 
8192(bottom). These have been boosted, downgraded to 
nside=512, and masked with a galactic mask provided by 
LAMBDA. Values are quoted in percent. 



FIG. 8. Comparison of nside=2048(top), 4096(middle), and 
8192(bottom). These have been boosted, downgraded to 
nside=1024, and masked with a galactic mask provided by 
LAMBDA. Values are quoted in percent. 



or 1024 give more accurate results than boosting and 
masking at the same resolution. 



VII. DISCUSSION 

We have investigated a real-space approach to deboost- 
ing a CMB temperature map by calculating fractional 
area overlaps of boosted pixels with the HEALPix grid in 
the spherical projection. We have shown that the resolu- 
tion at which the deboosting is performed greatly affects 



the resulting power spectrum and therefore one should 
take great care when extracting quantitative information 
from deboosted maps. We have shown that for a boost 
kernel that is calculated at Nside=8192, the highest res- 
olution used in this work, there is a 1% difference in the 
binned angular power spectrum at £ = 2000 between full- 
sky boosted and unboosted maps. Because boosting, ro- 
tating the map 180 degrees, and boosting again does not 
return the power spectrum to its unboosted value but in- 
stead further decreases the power at large multipoles at a 
2% level, the behavior is dominated by resolution effects 
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FIG. 9. Comparison of nside=2048(top) and 

nside=4096(bottom). These have been boosted and 
masked with a 20 degree isolatitude equatorial strip. Values 
are quoted in percent. 



and is not predominantly physical in origin. 

We have shown that resolution effects could lead one to 
incorrectly conclude that cosmological quantities were bi- 
ased by this systematic, as power spectra generated from 
CAMB with 1 sigma difference in n s and the WMAP best 
fit value of n s have percent differences at large multipoles 
consistent with those we see from boosted and unboosted 
maps at Nside=4096. 



Since masking maps mixes multipoles, we compared 
using masked maps that had boost contributions in- 
cluded to using unboosted masked maps for power spec- 
trum calculations. We found that, while masking along 
with boosting caused the difference between each indi- 
vidual Ci to vary more than for full-sky maps, masking 
boosted maps had little effect on the differences between 
binned Cg values. We found that if one wants to deboost 
a map and later mask with the publicly available galac- 
tic cuts released by WMAP (which can be downloaded 
at Nside=512 or Nside=1024) one should first deboost 
with a higher resolution boost kernel before masking the 
map. 

But, as Section |IV| shows, Lorentz boosts do move 
power to nearby multipoles - the behavior we see here is a 
feature of our particular power spectrum and the fact that 
it is a slowly varying function. It is not a generic state- 
ment about Lorentz boosts for any cosmological quan- 
tity and highlights a need for polarization fields and the 
three-point function to be investigated. 



ACKNOWLEDGMENTS 

We would like to thank Anthony Challinor, Arthur 
Kosowsky, Raul Abramo, Wayne Hu, Kendrick Smith, 
Maik Stuke, Dominik J. Schwarz, and Pascal Vau- 
drevange for useful conversations during the preparation 
of this work. TSP thanks the physics department of Case 
Western Reserve University and the Institute of Theoret- 
ical Astrophysics at Oslo for their hospitalities during the 
stages of this work. 

GDS and AY are supported by a grant from the US De- 
partment of Energy and were supported early in this work 
by NASA under cooperative agreement NNX07AG89G. 
AY is supported by NASA Headquarters under the 
NASA Earth and Space Science Fellowship Program - 
Grant TRN507323. 

We acknowledge the use of the Legacy Archive for Mi- 
crowave Background Data Analysis (LAMBDA). Sup- 
port for LAMBDA is provided by the NASA Office of 
Space Science. 



[1] M. Kamionkowski and L. Knox, "Aspects of the cosmic 
microwave background dipole," Phys. Rev. D 67, 063001 
(2003) [astro-ph/0210165| . 

[2] A. Challinor and F. van Leeuwen, "Peculiar velocity 
effects in high resolution microwave background ex- 
periments," Phys. Rev. D 65, 103001 (2002) [astro-| 
ph/01124"57] . 

[3] A. Kosowsky, T. Kahniashvili, "The Signature of Proper 
Motion in the Microwave Sky," |arXiv:1007.4539| [astro- 
ph.CO]]. 

[4] L. Amendola, R. Catena, I. Masina et at, "Measuring 
our peculiar velocity on the CMB with high-multipolc off- 
diagonal correlations," arXiv: 1008. 1183 [astro-ph.CO]]. 



[5] J. Chluba, "Aberrating the CMB sky: fast and accurate 
computation of the aberration kernel," arXiv: 1102.3415 
[astro-ph.CO]. 

[6] A. Notari and M . Quartin, JCAP 1202, 026 (2012) 
|arXiv:1112.l400l [astro-ph.CO]]. 

[7] R. Catena and A. Notari, "Cosmological parameter es- 
timation: impact of CMB aberration," arXiv: 1210.2731 
[astro-ph.CO]. 

[8] K. M. Gorski, E. Hivon, A. J. Banday, B. D. Wandelt, 
F. K. Hansen, M. Reinecke and M. Bartelman, Astro- 
phys. J. 622, 759 (2005) [astro-ph/0409513] . 

[9] McKinley, J.M., "Relativistic transformation of solid an- 
gle" , American Journal of Physics, 48, 612-614, 1980. 



[10] 



[11] 



T. S. Pereira, A. Yoho, M. Stuke et al, "Effects of a Cut, 
Lorentz-Boosted sky on the Angular Power Spectrum," 
jarXiv: 1009.49371 [astro-ph.CO]]. 

A. Lewis, A. Challinor, A. Lasenby, Anthony, "Effi- 
cient Computation of CMB anisotropics in closed FRW 
models", Astrophys. J., 538, 473-476, 2000. |astro-| 



ph/99lTT77l '. 

[12] G. Hinshaw et al. [WMAP Collaboration], "First year 
Wilkinson Microwave Anisotropy Probe (WMAP) obser- 
vations: The Angular power spectrum," Astrophys. J. 
Suppl. 148, 135 (2003) |astro-ph/03022T7]. 



